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Abstract. We present STAR results from identified particle spectra mea- 
sured in y/SNN = 62.4 GeV Au-Au collisions. Particle production and system 
dynamics are compared to results at ^snn = 200 GeV. We extract kinetic and 
chemical freeze-out parameters using blast wave model parameterization and 
', statistical model. We discuss the effect of resonance decays on the extracted 

f^) ■ kinetic freeze-out parameters. 
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■ Results from heavy ion collisions at RHIC generated great interest in the study 

of nuclear matter under extreme conditions. Theory predicts the existence of a 
new phase of matter at sufficiently high energy density which might have been 
reached in Au-Au collisions at RHIC. To test the hypothesis of a phase transition 
threshold and to map out the systematics, a lower energy Au-Au run has been 
performed at ^/snn = 62.4 GeV. In this paper we investigate the bulk particle 
production and dynamic evolution of the colliding system and contrast the 62.4 
GeV and 200 GeV results. We also address the kinetic and chemical freeze-out 
properties of the system and the effect of resonance decays on the extracted kinetic 
freeze-out parameters. 

The bulk of the produced particles (7r ± , A" ± , p and p) in heavy ion collision 
appears at low transverse momentum (px < 2 GeV/c). Bulk yield ratios of identified 
particles characterize the chemical freeze-out. Chemical freeze-out temperature, 
chemical potentials and the ad-hoc strangeness suppression factor are determined 
via statistical models. Kinetic freeze-out parameters - temperature and average 
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Fig. 1. Identified particle px spectra of negatively charged particles in 62.4 GeV 
Au-Au collisions. The lines indicate the blast wave model fits without including 
resonance decays. 



transverse flow velocity - are extracted from the shapes of identified particle spectra. 

2. Data analysis 

Measurements were carried out by the STAR experiment Charged particles are 
detected in the Time Projection Chamber (TPC) The data sample is divided 
into nine centrality classes according to the measured uncorrected charged particle 
multiplicity in \r)\ < 0.5. They represent 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 
40-50%, 50-60%, 60-70% and 70-80% of the total geometrical cross-section. 

Particles (tt^ , K ± , p and p) are identified by their specific energy loss (dE / dx) 
in the TPC gas. Spectra are corrected for acceptance, detector and tracking ineffi- 
ciencies, particle decays, and background processes. Corrections are obtained from 
Monte Carlo tracks embedded in real events at the raw data level with detailed 
description of the STAR detector geometry and realistic simulation of the TPC 
responses [0] El E]- Spectra presented here are inclusive measurements for X ± , p 
and p, while the pion spectra are corrected for weak decays. We followed the same 
analysis procedure as used in the measurement of identified particle spectra in 200 
GeV Au-Au collisions from STAR [0. 

3. Results 

Figure ^ shows the identified particle spectra of 7r ± , , p and p in the pr range 
of 0.2-1.2 GeV/c, dictated by the dE/dx capability. The evolution of the particle 
spectra is similar to that at 200 GeV [ E] : particle spectra harden with increasing 
centrality and with increasing particle mass because of collective transverse radial 
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Fig. 2. Left panel: Antiparticle-to-particle ratios at 62.4 GeV (solid markers) and 
200 GeV (empty markers) as a function of centrality. Right panel: K~ /tt~~ and 
p/jr~ ratios at 62.4 GeV (solid markers) and at 200 GeV (empty markers) as a 
function of centrality. 
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Fig. 3. Left panel: Extracted chemical (circles) and kinetic (triangles) freeze-out 
temperatures as a function of centrality. The 62.4 GeV data are shown in solid 
markers and 200 GeV data in empty markers. Right panel: Extracted average 
transverse flow velocities as a function of centrality. 



flow. 

Due to the lower collision energy, the net baryon density at mid-rapidity is 
increased as compared to 200 GeV. This is shown by the p/p and K~ /K + ratios in 
Fig. El (left panel). The p/p ratio decreases with centrality more rapidly at 62.4 GeV 
than at 200 GeV. Figure El (right panel) shows unlike particle ratios. Qualitatively 
similar centrality dependence is found at the two energies. Uncorrected charged 
multiplicity (Nch) is chosen here to indicate centrality, however, the qualitative 
features of the data are the same with other centrality measures, such as the collision 
geometry. 

Kinetic freeze-out parameters are extracted from the blast wave model without 
including resonances [|H|. The extracted kinetic freeze-out temperature {T^in) ex- 
hibits a decreasing trend with increasing centrality at both energies. The average 
transverse flow velocity ((/?)) increases with increasing centrality. These features 
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Fig. 4. Left panels: Calculated contributions to ir , K and p from resonance 
decays. Right panels: Resonance contributions relative to the primordial spectrum. 



are consistent with stronger expansion in more central collisions. 

The blast wave parameterization is used for extrapolation of the spectra to 
obtain the total yields. Particle yield ratios are used to extract chemical freeze- 
out parameters using local equilibrium models [ ED . The chemical freeze-out 
temperature (T c h) is independent of centrality and collision energy. The value of 
T c h is close to the QCD predicted critical temperature of phase transition [Ul lllllTS] . 

The extracted freeze-out parameters are summarized in Fig. [3] The similar 
results between 200 GeV and 62.4 GeV indicate similar collision evolution at both 
energies. 
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Fig. 5. Top left panel: Fit of the calculated spectra to the measured ones in top 
5% central Au-Au collisions at 200 GeV [\7\. Four calculated spectra are shown for 
7r~ (upper curves): including resonances with three different p contributions and 
excluding resonances. Only two calculated curves are shown for K~ (middle curves) 
and p (lower curves): including resonances with 100% p and excluding resonances. 
Other panels: data / calculation ratios. Error bars are from statistical and point- 
to-point systematic errors on the data, and are shown for only one set of the data 
points. 

4. Effect of resonance decays 

The kinetic freeze-out parameters above were extracted without including resonance 
decays (i.e. treating the measured particle spectra as primordial). Pion data below 
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375 MeV were excluded from the fit in order to minimize effect of resonance decays. 
In this section, we address the effect of resonance decays more rigorously. For 
illustration, the central (0-5%) 200 GeV Au-Au data is used. 

Our study is based on the model by Wiedemann and Heinz which cal- 

culates the pt spectra of thermal particles and decay products from resonances. 
There are, however, several distinctions in our study: (1) We assumed two distinct 
freeze-out temperatures while Ref. [ fTTTj used a single temperature. The relative 
proportions of particles and resonances in our study are determined by chemical 
freeze-out parameters that are fixed. We used T c h= 160 MeV, (ig = 22 MeV, Hs= 
1.4 MeV, and 7=0.98. In Ref. [EH the chemical potentials are set to zero. (2) We 
included a more extensive list of resonances, namely: p, u>, rj, rj', K*°, K^, <j> and 
A, A, S, 3, A1520, Si3 85 , f2. (3) Instead of Gaussian in Ref. [El, we used box 
profile for the flow velocity: (3 = j3s (r/R) n where n is fixed to be 0.82 [|7j. (4) 
We used a constant dN/dy instead of Gaussian in Ref. [ 13 . Rapidity distributions 
are needed for resonances which can decay into particles at mid-rapidity where our 
measurements are made. 

The resonance decay particle spectra are calculated, and combined with pri- 
mordial ones. Spin, isospin degeneracies and decay branching ratios are properly 
taken into account. 

The calculated particle spectra of tt~ , K~ and p are shown in Fig. ^ (left pan- 
els). Since the measured pion spectra are corrected for weak decays, the calculated 
inclusive pion spectra do not contain weak decay pions. Resonance contributions 
are labeled by the initial resonance particle, e.g. a tt emerging from the 77' — ^ 77 — )■ 7r 
decays is labeled as 7iy . The calculated inclusive pion spectra include contribu- 
tions from A1520 which are not plotted in Fig. 0] The right panels of Fig. 0] show 
the resonance contributions to the inclusive spectra relative to the primordial one. 
The low px 7t enhancement is the counter play of p, u> and 77; at higher px the p 
contribution dominates. The inclusive kaon and antiproton spectra do not show 
significant changes in the spectral shapes compared to the primordial ones. The 
largest contributions to the inclusive kaon spectra are from K*° and K*~ , and the 
largest contributions to the inclusive p and p spectra are from A, A, and S's. 

It is an open question what flow velocity and temperature should be assigned 
to the short lived resonances, such as p and A. These resonances are expected 
to be constantly regenerated during the system evolution, since their life-times are 
shorter than the expected system evolution time. 

In other words, the processes of p — ► ttit and tttt — > p, for example, constantly 
occur along the dynamical evolution of the system. Thus, it is reasonable to expect 
that the final p decay pions carry the same flow information as the primordial pions 
do. In other words, the regenerated p gain negligible flow velocity during its short 
life span except the inherited flow from the two resonant pions. 

To gain better insights, three cases are considered for p, which gives the largest 
contribution to the measured tt spectra. 

(1) The p decay pions have the same pr spectra shape as the primordial pions. 
(2) The p acquires flow as given by kinetic freeze-out temperature and transverse 
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flow velocity, and the decay pions are calculated from decay kinematics. (3) Half 
of the p contribution is taken like in (1) and the other half as in (2). Case (1) 
has the largest flow for decay pions because the p, being heavy, acquires flow more 
efficiently than pions. 

Tabled shows the fit results for the three cases. Also listed in Table Q] are the 
fit results without including resonances. The results are consistent with those in 
Rcf. [Ej. 

Figure shows the fits of the calculated inclusive spectra to the measured ones. 
Fits are performed to the six measured spectra simultaneously, but only negatively 
charged particles are shown. For K~ and p, results from the 100% p case fit and 
the fit excluding resonances are plotted, while all fits are plotted for pions. In case 
of 100% p, the calculated spectrum starts to deviate from data above pt ~ 400 
MeV. In case of 0% p, below pt ~ 400 MeV the calculated inclusive spectrum is 
enhanced by ui and 77, which become more important without the p. 

The model, with all the three cases of p contributions, seems to describe the data 
well. The spectra are found to be less sensitive to the kinetic freeze-out temperature 
than the flow velocity. Figure shows, as an example, the fitted x 2 versus fit 
parameters Tkin and ((3). It is seen from the figure that (/?) is well constrained 
but Tkin is less so. Nonetheless, the fitted Tki n values with all three cases of p 
contributions seem to agree with that obtained without including resonances within 
systematic error of ± 10 MeV [[?]. In other words, resonance decays appear to have 
no significant effect on the extracted kinetic freeze-out parameters as shown in 
Tabled This is primarily due to the limited pt ranges of our data where resonance 
decay products have more or less similar spectral shapes as the primordial particles 
do. 

We also fitted the spectra data with a single, fixed kinetic freeze-out tempera- 
ture Tkin = T c h = 160 MeV including resonances and with 100 % p contribution. 
The fitted {(3) is 0.520±o;oo2 witn X 2 /NDF=19.56. A single temperature scenario 
is therefore ruled out by the data. 

The x 2 /NDF is smaller than unity because we included in the fit the point-to- 
point systematic errors (dominate over statistical ones), which were estimated on 
the conservative side and might not be completely random. If we scale the x 2 /NDF 
such that the minimum is unity, then we get somewhat smaller statistical errors on 
the fit parameters. 

5. Conclusions 

In summary, we have reported results on identified particle spectra in 62.4 GeV 
Au-Au collisions and compared to those at 200 GeV. Significant drop is observed 
in the p/p ratio from 200 GeV to 62.4 GeV due to higher net baryon density at 
mid-rapidity at the lower energy. The extracted chemical freeze-out temperatures 
from statistical models show little dependence on collision centrality or energy, and 
coincide with the predicted phase transition temperature. The extracted kinetic 
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Fig. 6. Fitted x 2 as a function of Tkm and (/?) for the 0% p case. Number of 
degrees of freedom NDF — 75. 



freeze-out temperature drops with centrality while the radial flow velocity increases 
with centrality, consistent with stronger expansion in more central collisions. The 
freeze-out properties show similar evolution at both energies. 

We carried out a rigorous study of resonance decay effect on the extracted 
kinetic freeze-out parameters from the blast wave model. With different contribu- 
tions from the short lived p resonance, the extracted parameters agree with those 
obtained without including resonance decays within systematic uncertainties. Res- 
Table 1. Extracted kinetic freeze-out parameters and fit \ 2 ■ The flow profile n 
parameter is fixed to 0.82. 

~Set T km (MeV) (/3) X 2 /ndf 



No resonances 86.8_ ; 6 

p % 94.6±H 

p 50 % 87.41?;? 

p 100 % 77.2±g;| 



0.595™ 0.26 

0.603™ 0.37 

0-605™ 0.45 

0.60T™ o.60 
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onance decays appear to have no significant effect on the extracted parameters. 
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